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Abstract. The Genetic Algorithm is a heuristic that can be used to produce model 
independent solutions to an optimization problem, thus making it ideal for use in cosmology 
and more specifically in the analysis of type la supernovae data. In this work we use the Genetic 
Algorithms (GA) in order to derive a null test on the spatially flat cosmological constant model 
ACDM. This is done in two steps: first, we apply the GA to the Constitution SNIa data in 
order to acquire a model independent reconstruction of the expansion history of the Universe 
H(z) and second, we use the reconstructed H(z) in conjunction with the Om statistic, which 
is constant only for the ACDM model, to derive our constraints. We find that while ACDM is 
consistent with the data at the 2a level, some deviations from ACDM model at low redshifts 
can be accommodated. 



1. Introduction 

In the previous decade it was discovered that the Universe is undergoing an accelerated expansion 
[I] j [2] j [2] ■ This acceleration is usually attributed either to a cosmic fluid with negative pressure 
dubbed Dark Energy or to an IR modification of gravity. In order to identify the properties of 
Dark Energy or the structure of the IR modification of gravity it is necessary to know to a high 
precision the rate of the expansion of the Universe, parameterized as H = ~ where a = j^— is the 
scale factor and z is the redshift of the cosmological probe, as it measured by the observations. 

The behavior of the expansion of the Universe can be identified by studying two functions, 
the Equation of State (EoS) w(z) = which can be rewritten as 

*(*) = -1 + 3(1+*) v , (1) 

where 5H(z) 2 = H(z) 2 /Hq — f2om(l + z ) 3 accounts for all terms in the Friedmann equation not 
related to matter and the deceleration parameter q(z) = which can be rewritten as 

, W = -l + (l + 2 )^M§(f», (2) 

Obviously, the cosmological constant (w(z) = — 1) corresponds to a constant dark energy density, 
while in general w{z) can be time dependent. Also, an important parameter is the value of the 
deceleration parameter today, ie q(z = 0) = qo, which for the cosmological constant model in 
GR it is q = -1 + 3ft m/ 2 - 

However, despite all the recent progress the origin of the accelerated expansion of the universe 
still remains unknown with many possibilities still remaining open, see for example [3]. The 



simplest choice that agrees well with the data is a positive cosmological constant which has to 
be small enough to have started dominating the universe at late times. As it was demonstrated 
by the Seven- Year WMAP data [5], the cosmological constant remains the best candidate and 
has the advantage of having only one free parameter related to the properties of the Dark Energy. 
Nonetheless, this model fails to explain why the cosmological constant is so small that it can 
only dominate the universe at late times, a problem known as the coincidence problem and there 
are a few cosmological observations which differ from its predictions [6J,[7J. 

A very important complication in the investigation of the behavior of dark energy occurs due 
to the bias introduced by the parameterizations used. At the moment, there is a multitude of 
available phenomenological ansatze for the dark energy equation of state parameter w or dark 
energy density, each with its own merits and limitations (see [8J and references therein). The 
interpretation of the SNIa data has been shown to depend greatly on the type of parametrization 
used to perform a data fit [SJ , [9J . Choosing a priori a model for dark energy can thus adversely 
affect the validity of the fitting method and lead to compromised or misleading results. 

The need to counteract this problem paved the way for the consideration of a complementary 
set of non-parametric reconstruction techniques [ID] . [IT], [12], [13], [H] and model independent 
approaches ([E],[l6], [9J, [IT], [IS], [15]. These try to minimize the ambiguity due to a possibly 
biased assumption for w by fitting the original datasets without using any parameters related 
to some specific model. The result of these methods can then be interpreted in the context of a 
dark energy model of choice. Non-parametric reconstructions can thus corroborate parametric 
methods and provide more credibility. However, they too suffer from a different set of problems, 
mainly the need to resort to differentiation of noisy data, which can itself introduce great errors. 

In this paper, we present a method that can be used as a model independent approach in 
testing the standard cosmological model. This is done by using the Genetic Algorithms (GA) 
technique, first used in the analysis of SNIa data in [20]. The GAs represent a method for 
non-parametric reconstruction of the dark energy equation of state parameter w, based on the 
notions of genetic algorithms and grammatical evolution. GAs are more useful and efficient than 
usual techniques when 

• The parameter space is very large, too complex or not enough understood, as is the case 
with dark energy. 

• Domain knowledge is scarce or expert knowledge is difficult to encode to narrow the search 
space. 

• Traditional search methods give poor results or completely fail. 

Naturally, therefore, they have been used with success in many fields where one of the above 
situations is encountered, like the computational science, engineering and economics. Recently, 
they have also been applied to study high energy physics [21], [22], [23] gravitational wave 
detection [24] and gravitational lensing [25J. Since the nature of Dark Energy still remains 
a mystery, this makes it for us an ideal candidate to use the GAs as a means to analyze the 
SNIa data and extract model independent constraints on the behavior of the Dark Energy. 

This talk was delivered at the 14th conference in the series "Recent Developments in Gravity" 
(NEB- 14) and is a brief overview of the use of the Genetic Algorithms to a model independent 
reconstruction of the expansion history of the Universe. In Section 2 a brief description of the 
Om statistic is given while in Section 3 we provide an overview of the general methodology of 
the GA paradigm and finally, we present our results in Section 4. 

2. The Om statistic 

The recently introduced Om diagnostic [26J (also look at Ref. [27]) enables us to distinguish 
ACDM from other dark energy models without directly involving the cosmic EOS. The Om 
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Figure 1. The Om diagnostic is shown as a function of redshift for dark energy models with 
n 0m = 0.27 and w = —1,-0.8,-0.6,-0.4,-0.2 (bottom to top). For Phantom models (not 
shown) Om would have the opposite curvature. 



diagnostic is defined as: 



Om(x) = 



h 2 (x) - 1 
3-1 ■ 
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l + z 
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h(x) = H{x)/Hq . 



For dark energy with a constant equation of state w = const, 

h 2 (x) = n 0m x 3 + (1 - n 0m )x a , a = 3(1 + w) 
(we assume that the universe is spatially flat for simplicity). Consequently, 

,x a - 1 



Om{x) = Qom + (1 — Jl 



0m) 



1 ' 



(3) 



(4) 



(5) 



from where we find Om(x) = $7om in ACDM, whereas Om(x) > Q,Q m in quintessence (a > 0) 
while Om(x) < Qo m in phantom (a < 0). We therefore conclude that: Om(x) — Ho m = if dark 
energy is a cosmological constant. Note that since S1a + Slom — 1 in ACDM, this model contains 
SCDM (r^om = 1,^A = 0) as an important limiting case. Consequently, the Om diagnostic 
cannot distinguish between large and small values of the cosmological constant unless the value 
of the matter density is independently known. In other words, the Om diagnostic provides us 
with a null test of the cosmological constant. This is a simple consequence of the fact that h 2 (x) 
plotted against x 3 results in a straight line for ACDM, whose slope is given by fiom- For other 
dark energy models the line describing Om(x) is curved, since the equality 



d[h 2 ] 
d[x 3 } 



constant , 



(6) 



(which always holds for ACDM for any x = 1 + z) is satisfied in quintessence/phantom type 
models only at redshifts significantly greater than unity, when the effects of dark energy on the 
expansion rate can safely be ignored. As a result the efficiency of the Om diagnostic improves 
at low redshifts (z < 2) precisely where there is likely to be an abundance of cosmological data 
in the coming years! 

For a constant EoS: 1 + w ~ [Om(z) — Slom](l — ^Om) -1 at z <C 1, consequently a larger 
Om(z) is indicative of a larger w; while at high z, Om(z) — > £lom: a s shown in Fig. [TJ 

On the other hand, for quintessence as well as phantom the line describing Om(x) is curved, 
which helps distinguish these models from ACDM even if the value of the matter density is not 
accurately known - see Fig. [TJ 

In practice, the construction of Om requires a knowledge of the Hubble parameter, h(z), 
which can be determined using a number of model independent approaches [IS]) [IS]) [2]) [T7] . 

mm- 

3. Genetic algorithms 

3.1. Overview 

GAs were introduced as a computational analogy of adaptive systems. They are modeled loosely 
on the principles of the evolution via natural selection, employing a population of individuals that 
undergo selection in the presence of variation-inducing operators such as mutation and crossover. 
The encoding of the chromosomes is called the genome or genotype and is composed, in loose 
correspondence to actual DNA genomes, by a series of representative "genes". Depending on 
the problem, the genome can be a series (vector) of binary numbers, decimal integers, machine 
precision integers or reals or more complex data structures. 

In order to evaluate the individual chromosomes a fitness function is used and reproductive 
success usually varies with fitness. The fitness function is a map between the gene sequence of 
the chromosomes (genotype) and a number of attributes (phenotype), directly related to the 
properties of a wanted solution. Often the fitness function is used to determine a "distance" of 
a candidate solution from the true one. Various distance measures can be used for this purpose 
(Euclidean distance, Manhattan, Mahalanobis etc.). 

The algorithm begins with an initial population of candidate solutions, which is usually 
randomly generated. Although GA's are relatively insensitive to initial conditions, i.e. the 
population we start from is not very significant, but using some prescription for producing this 
seed generation can affect the speed of convergence. In each successive step, the fitness functions 
for the chromosomes of the population are evaluated and a number of genetic operators (mutation 
and crossover) are applied to produce the next generation. This process continues until some 
termination criteria is reached, e.g. obtain a solution with fitness greater than some predefined 
threshold or reach a maximum number of generations. The later is imposed as a condition to 
ensure that the algorithm terminates even if we cannot get the desired level of fitness. 

The various steps of the algorithm can be summarized as follows: 

(i) Randomly generate an initial population M(0) 

(ii) Compute and save the fitness for each individual m in the current population M(t). 

(iii) Define selection probabilities p(m) for each individual m in M(t) so that p(m) is proportional 
to the fitness. 

(iv) Generate M(t+ 1) by probabilistically selecting individuals from M{t) to produce offspring 
via genetic operators (crossover and mutation). 

(v) Repeat step 2 until a satisfying solution is obtained, or a maximum number of generations 
reached. 

We should stress that the initial population M(0) will only depend on the choice of the 
grammar and therefore it can only affect on how fast the algorithm will converge to the minimum. 



Using the wrong grammar may result in the GA being trapped in a local minimum. Also, two 
important parameters that affect the transition from the population M(t) to M(t + 1) are the 
selection rate and the mutation rate. The selection rate is typically of the order of 10% and it 
affects how many of the individuals, after they are ranked with respect to their fitness, will be 
allowed to produce offspring. The mutation rate is usually of the order of 5% and it expresses 
the probability that an arbitrary part of the genetic sequence will be changed from its previous 
state. This is done in order to maintain the genetic diversity from one generation of a population 
to the next. 

The paradigm of GAs described above is usually the one applied to solving most of the 
problems presented to GAs. Though it might not find the best solution, more often than not, it 
would come up with a partially optimal solution. A more detailed overview and its application 
in cosmology can be found in |2Uj . 

The difference of the GA over the the standard analysis of the data, ie having an a priori 
defined model with some free parameters, is that the obtained values of the best-fit parameters 
are model-dependent and in general models with more parameters tend to give better fits to the 
data. This is where the GA approach starts to depart from the ordinary parametric method. Our 
goal is to minimize a function, not using a candidate model function for the distance modulus 
and varying parameters, but through a stochastic process based on a GA evolution. This way, 
no prior knowledge of a dark energy model is needed to obtain a solution and our result will be 
completely parameter- free. 

In simpler terms, the GA method does not require an a priori assumption for a DE model, 
but uses the data themselves to find this model. Also, it is parameter free as the end result does 
not have any free parameters that can be changed in order to fit the data. So, in this sense this 
method has far less bias any of the standard methods for the reconstruction of the expansion 
history of the Universe. This is the main reason for the use of the GAs in this paper. 

3.2. General Methodology 

We first outline the course of action we follow to apply the GA paradigm in the case of SNIa 
data. For the application of GA and grammatical evolution (GE) on the dataset, we use a 
modified version of the GDF [25] tool 0, which uses GE as a method to fit datasets of arbitrary 
size and dimensionality. This program uses the tournament selection method for crossover. GDF 
requires a set of input data (train set), an (optional) test sample and a grammar to be used 
for the generation of functional expressions. The output is an expression for the function which 
best fits the train set data. 

Since the SNIa datapoints are given in terms of the distance modulus /x ^ s (zi), the fitness 
function for the GA we have chosen is equal to — XsNia> wnere 

2 _ (Vobs(zj) - LLGA(Zj)) 2 , , 

XsNIa — 72 ' \') 

i=l V i 

where for the Constitution set N = 397 and a 2 i are the errors due to flux uncertainties, intrinsic 
dispersion of SNIa absolute magnitude and peculiar velocity dispersion and HGA(zi) is the reduced 
distance modulus obtained for each chromosome of the population by the GA. 

An advantage of our method is the fact that neither the XsNIa nor ^ e Hga(z) depend on 
any parameters. Also, as it can be seen from eq. (|7j), our method does not make an assumption 
of flatness or not in order to fit the data. Flatness (or non-flatness) only comes into play when 
one tries to find the luminosity distance di,(z) and the underlying dark energy model given the 
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best-fit form of /jiga(z)- So, in this aspect the result produced by the GA is independent of the 
assumption of flatness as well. 

The GA evaluates © in each evolutionary step for every chromosome of the population. 
The one with the best fitness will consequently have the smallest XsNia an d De the best 
candidate solution in its generation. Of all the steps in the execution of the algorithm, the 
evaluation of the fitness is the most expensive. GDF is appropriately modified to use ([7]) as the 
basis of its fitness calculation. 

After the execution of the GA, we obtain an expression ^lga{z) for the reduced distance 
modulus as the solution of best fitness and the corresponding XsNIa- Using this we can, through 
differentiation, obtain other functions or parameters of interest, such as the Om statistic or the 
deceleration parameter q(z). For example, for a flat Universe the Hubble parameter will be 
given by 

and then the deceleration parameter can be found by Eq. (|2J). 

A possible complication of this method is that it gives no direct way to estimate the errors 
for the derived parameters. One cannot expect to get an estimate by just running the algorithm 
many times and obtaining slightly different parameters. The GA usually tends to converge at 
the same solution for a given dataset, unless we change significantly the population size or the 
number of generations. A way to circumvent this problem is to use a bootstrap Monte Carlo 
simulation to produce synthetic datasets and rerun the algorithm on them. We can thus obtain 
a statistical sample of parameter values which will allow us to estimate the error. 

We sketch the procedure we will follow 



(i) The GA is applied on the original SNIa dataset with the chosen execution parameters and 
a solution for P>ga(z) is obtained. 

(ii) Generate a number of synthetic datasets by drawing each time the same number of data 
points with replacement from the original set. 

(iii) The GA is rerun for the synthetic datasets. 

(iv) A new set of values for the Om statistic and the deceleration parameter q(z) is generated. 

(v) The 95% error of the desired parameter can be found by taking the 2.5 and the 97.5 
percentiles of the bootstrap distribution [30] . 

Using the above steps, we can obtain error estimates for any desired parameter. 



3.3. An example 

In this Section we will briefly describe a simple exampk^] of how he GA determines the best-fit. 
We will avoid describing the technicalities, like the binary representation of the solutions, and 
instead we will concentrate on how the generations evolve, how many and which individuals are 
chosen for the next generation etc. 

As we mentioned earlier, the choice of the grammar is very important, however in order to 
keep our example as simple as possible we will assume that our grammar includes only basic 
functions like polynomials x, x 2 etc, the trigonometric functions sin(x), cos(x), the exponential 
e x and the logarithm ln(x). 

The first step in the GA is setting up a random initial population M(0) which can be 
any simple combination of these functions, eg Hga,i{z) = ln(z), /j*ga,2(z) = —1 + z + z 2 and 



2 This is only a schematic description of how the GA works and this is done solely for the sake of explaining the 
basic mechanisms of the GA. 



fiGA 3(2) = sin(z). The number of number of candidate solutions (chromosomes) in the genetic 
population is usually a few hundreds and later on we will use the value of 500. 

Next, the algorithm measures the fitness of each solution by calculating their x 2 \ i e f° r 
our simple example xl = 80579.9, xl = 292767.0 and xl = 412928.0. The selection per se 
is done by implementing the "Tournament selection" method which involves running several 
"tournaments" , sorting the population with respect to the fitness each individual and after that 
a fixed percentage, adjusted by the selection rate (see the previous subsection), of the population 
is chosen. As we mentioned, the selection rate is of the order of 10% of the total population, 
but lets assume that the two out of the three candidate solutions QucM,i (z) and Hga,2{z)) are 
chosen for the sake of simplicity. 

The reproduction of these two solutions will be done by the crossover and mutation operations. 
The crossover will randomly combine parts of the "parent" solutions, for example in one such 
realization this may be schematically shown as 

VGA,l(z) ®VGA,2{z) -> (fiGA,l(z),P>GA,2(z),fiGA,3( z )) 

= (ln(z 2 ),-l+ln(z 2 ),-l + ln(z)) 

After this is done, the GA will proceed to implement the mutation operation. The probability 
for mutation is adjusted by the mutation rate, which as we mentioned is typically of the order 
of 5%. In our example this may be a change in the power of some term or the change in the 
value of a coefficient. For example, for the candidate solution Jjlgaz{ z ) this can be schematically 
shown as p,GA,3( z ) = — 1 + ln(z) — >• — 1 + ln(z 3 ), where the power of the z term was mutated 
from 1 to 3. 

Finally, at the end of the first round we have the three candidate solutions 

M ( l ) = (pGA,l{z),p,GA,2(z), fiGA,3(z)) = 
(ln(z 2 ), -1 + ln(z 2 ), -1 + ln(z 3 )) 

At the beginning of the next round the fitness of each individual will again be determined, 
which for our population will be (x 2 , xl,xl) = (3909.2,18435,113327), and the selection and 
the other operations will proceed as before. It is easily seen that even after one generation the 
X 2 of the candidate solutions has decreased dramatically and as we will see in the next section, 
it only takes a bit more than 100 generations to reach an acceptable x 2 and about 500 to 1000 
generations to start converging on the minimum. 

After a predetermined number of generations, usually on the order of 10 3 , has been reached 
or some other criteria have been met, the GA will finish, with the best candidate solution of the 
last generation being the best-fit to the data. 

4. Results 

We will now examine a number of different configurations, highlighted in Table [TJ In order to 
keep our analysis as simple as possible we will use a Basic Grammar that includes polynomials 
like x, x 2 etc, the trigonometric functions sin(x), cos(x), the exponential e x and the logarithm 
ln(x) and a second grammar that includes the Basic Grammar and the Legendre polynomials. 
For the analysis we will use the Constitution SNIa dataset of [31J consisting of 397 SNIa. The 
steps followed for the usual minimization of ([7]) in terms of its parameters are described in detail 
in [32], [33], [S], [35]. 

In Fig. [2] we show the fitness (equal to — x 2 ) as a function of the number of generations. 
Clearly, all configurations seem to have to have converged very early or to be very close to 
converging to their minimum. An exception is Case 3 which even after 3000 generations still 
has not converged and we had to run the simulation for another 3000 generations to check its 



Table 1. The different cases considered in the analysis. The Basic Grammar includes 
polynomials, trigonometric functions, exponentials and logarithms. For ACDM the best-fit 
corresponds to VL m = 0.289. 



Case 1 


Basic 


grammar 


y 2 • 


= 465.19 


Case 2 


Basic grammar 


+ Legendre polyn. 


A.rmra 


= 477.87 


Case 3 


Only polyn. in £ 


grammar (3000 gen.) 


y 2 • 


= 468.19 


Case 4 


Only polyn. in £ 


grammar (6000 gen.) 


y 2 • 


= 462.56 


ACDM 


h(z) 2 = H$(n m (i + z) 3 + i- n m ) 


y 2 ■ 


= 465.51 



convergence. Since we want to compare the configurations with equal standards, ie number of 
generations etc, we will treat this case (Case 4) separately. 

As it can be seen from Table Q] and Fig. [31 most configurations agree well with the data and 
actually have a minimum y 2 comparable or even better with respect to ACDM. An exception 
to this is Case 2, where the best fit differs from ACDM by a Ay 2 = 12 thus providing a much 
poorer fit. The reason for this is the fact that by increasing the grammar we effectively increase 
the landscape in which the algorithm has to search and therefore, we are making the fit poorer 
as the GA will take much much longer to converge. 
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Figure 2. The fitness (equal to — y 2 ) 
function of the number of generations. The 
magenta line corresponds to case 1, the green 
line to case 2 and the blue line to case 3. The 
inset graphic shows case 4. 



Figure 3. The SNIa distance modulus against 
the redshift. The red line represents the best- 
fit ACDM model with Vt m = 0.289 while the 
magenta, blue, green and cyan (best-fit) lines 
correspond to cases 1-4 respectively. 



In Figs. H] and [5] we show the corresponding results of each case for the Om statistic and 
the deceleration parameter q{z) (our best-fit function corresponds to the cyan line). Due to 
the presence of trigonometric functions into the best fit and the fact that the Om statistic is 
derived from the Hubble parameter by differentiation, Case 1 exhibits strong oscillations. This 
is also apparent in the corresponding plot for the deceleration parameter for Case 2 (magenta 
line, Fig. [5]). 

As it can be seen in Fig. [H a constant w (but different from w = —1) forces the Om statistic 
to be a monotonic function that is located completely either above or below the constant case 



(Om = 0.289). However, none of the curves in Fig. H]are monotonia functions and all of them 
cross the constant Om case several times. Also, all cases remain reasonably close to the ACDM 
value of tt m = 0.289 for intermediate redshifts but deviate from that value at high or low 
redshifts. Therefore, this could be the effect of a time-dependent equation of state w(z) instead 
of a constant w. Comparing with Fig. [TJ it is possible that this behaviour is due to an equation 
of state w(z) more negative than -1. 

This is also supported from Fig. H] and Fig. as at low redshifts all cases predict a much 
different forms of acceleration than the best fit ACDM. In Fig. [3] and Fig. \5\ our best fit results 
show an overall slowing down of the acceleration at the low redshifts (which is in concordance 
with results presented in [S] with a very rapid increase of acceleration at the very low redshifts. 
These are very similar with the results presented in Fig. 2 of Ref. [TS] where a smoothing method 
for the reconstruction of the expansion history of the universe is being used. It is obvious that 
when two very different model independent approaches reach to a similar result, we must be on 
the right track and very close to the best possible result. However, theoretical interpretation of 
these results is beyond the scope of this paper. In fact there might be very different mechanisms 
that result to a similar expansion history of the universe as we reconstructed in this paper (e.g 
see [36], [E], [37], [38] ). 

In Fig. [6] we show the Om statistic as a function of the redshift z. The black line corresponds 
to the best fit of case 4, while the gray-shaded area to the 2a error region. The error region 
was calculated by implementing a bootstrap monte-carlo simulation as discussed in the previous 
section. Note that the overall slowing down of the acceleration at the low redshifts cannot be 
clearly seen from this figure as the errors are quite symmetrical around the best fit ACDM 
value. However, as we mentioned earlier this is clearly seen by the behavior of the deceleration 
parameter in Fig. [5j 

As it can be seen in Fig. [6l ACDM remains consistent with the data at the 2a level but at 
the same time, some deviations from the standard ACDM model can be accommodated. Due 
to the large errorbars of the current datasets, many dark energy models are still consistent with 
them but we expect this to change in the near future with the advent of high quality Snla 
observations. 

In Fig. [7J we show a histogram of the bootstrap distribution found from the monte-carlo 
simulation that was used to create the error regions. Clearly, the distribution is quite 
symmetrical and centered, so we are confident that our error region in Fig. [6] was correctly 
reconstructed. Finally, we also tested the GA on a mock ACDM dataset and we found that it 
can reconstruct the luminosity distance curve to within 0.1%, so we are confident that is does 
not overfit the data, but instead it can successfully detect the underlying DE model. 

5. Conclusions 

We used the Om statistic and the GAs in order to derive a null test on the cosmological constant 
model ACDM. Our interest in the GAs stems from the fact that they represent a method for 
non-parametric reconstruction of the expansion history of the universe, based on the notions 
of genetic algorithms and grammatical evolution. These kinds of algorithms are more useful 
and efficient than usual techniques especially when the problem under examination is not well 
understood, as is the case with dark energy. Since the nature of dark energy still remains a 
mystery, this makes it for us an ideal candidate to use the GAs as a means to analyze the SNIa 
data and extract model independent constraints on the behavior of the Dark Energy. On the 
other hand, the Om diagnostic (Sahni et al. 2008) enables us to distinguish ACDM from other 
dark energy models without directly involving the cosmic EoS. 

Our methodology is completely model independent and it can be summarized in two steps: 
first, we applied the GA to the Constitution SNIa data in order to acquire a model independent 
reconstruction of the expansion history of the Universe H(z). After we reconstructed H(z), the 




Figure 4. The Om statistic against the redshift. 
The black dashed line represents the ACDM 
value of f2 m = 0.289 while the magenta, blue, 
green and cyan (best-fit) lines correspond to 
cases 1-4 respectively. 



Figure 5. The deceleration parameter q against 
the redshift. The red line represents the best- 
fit ACDM model with VL m = 0.289 while the 
magenta, blue, green and cyan (best-fit) lines 
correspond to cases 1-4 respectively. 
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Figure 6. The Om statistic as a function of z. 
The black line corresponds to case 4 (the best- 
fit) and the gray-shaded area to the 2a error. 



Figure 7. Histogram of the bootstrap distribu- 
tion found from the monte-carlo simulation that 
was used to create the error regions of Fig. [6) 



second step was to use it in conjunction with the Om statistic and derive our null test. 

Our main result was that ACDM remains consistent with the data at the 2a level, see Fig. [6j 
but at the same time, some deviations from the standard ACDM model can be accommodated 
and this is also in accordance with the data especially at low redshifts (see the behavior of the 
best fit in Fig. [6]) . 

However, it should be mentioned that the slowing down at low redshifts for the Constitution 
set mentioned earlier can also be seen with the standard analysis with the CPL ansatz (see |39j) 
but only for the Constitution set. For the other SNIa datasets the reverse trend (speeding up 
the acceleration at low z) was observed. For this reason and due to the fact that the current 
data still have quite large errors we believe that the full potential of our method, while it is 
very promising, will only be realized in the near future when more high quality SNIa data will 
become available. When this happens, the error region at low z in Fig. [6] will be small enough 
for the Om statistic to completely discriminate between ACDM and the various Dark Energy 



models. 
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